Two-band model of Raman scattering on iron pnictide superconductors 
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Based on a two-band model, we study the electronic Raman scattering intensity in both normal 
and superconducting states of iron-pnictide superconductors. For the normal state, due to the match 
or mismatch of the symmetries between band hybridization and Raman vertex, it is predicted that 
overall B lg Raman intensity should be much weaker than that of the Bi g channel. Moreover, in the 
non-resonant regime, there should exhibit a interband excitation peak at frequency uj ~ 7.3ii(6.8£i) 
in the P>i g (B29) channel. For the superconducting state, it is shown that /3-band contributes most 
to the B2 g Raman intensity as a result of multiple effects of Raman vertex, gap symmetry, and 
Fermi surface topology. Both extended s- and d xy -wa.ve pairings in the unfolded BZ can give a good 
description to the reported P>2 S Raman data [Muschler et ai, Phys. Rev. B. 80, 180510 (2009).], 
while d x 2_ y 2-wave pairing in the unfolded BZ seems to be ruled out. 

PACS numbers: 74.20.-z, 74.20.Gz, 74.25. Jb 
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I. INTRODUCTION 

Recently superconductivity has been observed in sev- 
eral classes of iron-pnictide materials — ~— Despite that dif- 
ferent classes of pnictides have somewhat different crystal 
structures, their electronic structures behave to be quite 
similar, as confirmed by ARPES and magneto-oscillation 
measurements!*— Band structure calculations indicate 
that iron pnictides exhibit a quasi-two-dimensional elec- 
tronic structure. Their parent compounds are metals 
which display an antiferromagnctic long-range order. Su- 
perconductivity is induced cither by hole doping or elec- 
tron doping when part of Fe 2+ ions are replaced by Fe + , 
or solely by the application of high pressure. Overall 
speaking, iron-pnictide shows a good resemblance with 
the high-T c cuprate superconductor, and hence is an ideal 
candidate for studying the superconducting (SC) mech- 
anism of high T c . 

One crucial issue towards identifying the interaction 
that drives superconductivity is the symmetry structure 
of Cooper pairs. A conclusive observation of the pair- 
ing symmetry remains unsettled for iron pnictides, how- 
ever, Both nodal and nodclcss order parameters were re- 
ported in experimental observations. ARPES measure- 
ments clearly indicated a nodclcss gap at all points of the 
Fermi surface— ^ Moreover, magnetic penetration depth 
measurement has revealed a T 2 behavior down to 0.02T C 
- a possible signature for the unconventional s± state— 
The s± state is currently a promising pairing candidate 
for iron pnictides, which has a sign reversal between 
a and ft bands and can be naturally explained by the 
spin fluctuation mechanism—"— On the other hand the 
neutron-scattering spectrum seems to suggest the order 
parameter being fully gapped but without a sign rever- 
sal (called s++-wave in contrast to s-j--wave)i 12 i 13 On the 
contrary, both scanning SQUID microscopy^ and NMR 
measurement 1 5 d 6 are in favor of nodal SC order parame- 
ters. For phase sensitive experiments however, one point- 
contact spectroscopy reported was in favor of a nodal 



gapi 7 -, while the others reported were in favor of a node- 
less gap— - 

Compared to other spectroscopy experiments, Raman 
scattering has a unique feature, namely, by manipulat- 
ing the polarization directions of incident and scattered 
photons with respect to the crystallographic directions, 
it can selectively excite quasiparticles (QPs) on differ- 
ent parts of the Brillouin zone (BZ) . Therefore Raman 
scattering is considered a very useful probe for study- 
ing pairing symmetry of the SC order parameter. More 
explicitly, among various symmetry channels (A\ g , A2 g , 
B\g, and Bi g within a symmetry group), powers 
laws of low-energy Raman intensity as well as peak po- 
sitions can often give useful information on the pair- 
ing symmetry of a superconductor. Raman scattering 
had been very successful on the studies of the high-Tc 
cuprate superconductors—^ It is also anticipated that Ra- 
man scattering will give very useful signals on the iron- 
pnictide superconductors. 

Electronic Raman scattering experiment has re- 
cently been carried out on high-quality single- 
crystal Ba(Fei_ a; Co a ;)2As2.— Both normal- and 
superconducting-state data were taken and for fre- 
quency below 300 cm -1 , it showed that a significant 
change of the Raman response between the normal and 
superconducting state occurs only in the Bi g channel 
to which a strong peak develops at lu = 69cm _1 in the 
superconducting state. Moreover, the low-energy spec- 
tral weight in the SC state suggests a gap having nodes 
for certain doping levels of Ba(Fci_ x Co x )2As2. It is of 
importance to perform a more quantitative theoretical 
fitting on the Raman data reported in Ref. [2(| and this 
is indeed the goal of the present paper. 

On the theoretical side, Boyd et al£L have calculated 
the Raman response for iron-pnictides taking into ac- 
count multiple gaps on different Fermi sheets. In their 
calculations, the Raman vertices are obtained by the ex- 
pansion of harmonic functions for a cylindrical FS and 
the screening effect due to the long-range Coulomb in- 
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teraction is considered but leaving out vertex correc- 
tions capturing effects of final state interactions. Their 
results give a criterion to distinguish the momentum 
and frequency dependence of the SC order parameters 
by Raman scattering. In another study, Chubukov and 
coworker a 22 ' 23 have also calculated the Raman response 
for iron-pnictide superconductors. By analyzing the ver- 
tex corrections for the extended s-wave gap symmetry, 
they have revealed a collective mode below 2Aq in the 



wavevector k, respectively. Energy dispersions of orbital 
d xz and d yz are given by e x = —2t\ cos k x — 2i 2 cos k y — 
4f3 cos k x cos k y and e y = —2t2COsk x — 2ticos£; y — 
4^3 cos k x cos k y respectively, which are coupled by the 
d xy orbital with dispersion e xy = — 4t^ sin k x sin k y . By 
the unitary transformation 



cos 6>k — sm #k 
sin (9k cos 6*k 



ak<7 

Pka 



A lg channel that may be crucial to unravel the pairing Hamiltonian p can be transformed to 



symmetry. It provides an alternative way to distinguish 
between various suggested gap symmetries of the iron- 
pnictide superconductors. 

The electronic structure of iron-pnictide materials is 
considered to be more complex than that of the high- 
T c cuprates. The band structure calculations showed 
that superconductivity is associated with the Fe-pnictide 
layer, and the FS consists of two hole pockets and two 
electron pocketsi 24 ' 25 The maximum contribution to the 
density of states near the FS is due to the Fe-3d or- 
bitals. Several tight-binding models are thus proposed 
to construct the band structure and several competing 
orders are suggested to exist in these materials i^Sr— In 
this paper, we shall study the Raman response based on 
a two-band model.— This minimal model is considered to 
be a promising one which captures the major features of 
the electronic structures of several classes of iron-pnictide 
materials. As mentioned before, our aim is to give a more 
quantitative fitting to the current available Raman data 
and uncover useful fitting parameters. It is hoped that 
our st udy , together with other theoretical Raman works 
(Refs. l2ll andl22l) can give a more complete picture on the 
pairing symmetry and electronic structure of iron pnic- 
tides. 

The paper is organized as the following. In Sec. [TT1 
we present the two-band model Hamiltonian for iron- 
pnictide superconductors. In Sec. IIII| theoretical deriva- 
tions are given for the Raman scattering both in the nor- 
mal and in the SC states with respect to the two-band 
model given in Sec. HU Calculations of Raman response 
are presented for both normal and SC states. Some pre- 
dictions are made. Especially fitting to the Bi g Raman 
peak in the SC state are performed with useful fitting 
parameters given. Sec. IIVI is a summary of this paper. 



II. MODEL HAMILTONIAN 

For iron-pnictides, we start from the minimal two-band 
model introduced in Ref. [29[. The normal-state Hamil- 
tonian reads as 
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where E± = s+ ± J e 2 _ + e 2 y — /i with e± = (e x ± e y ) /2. 

The new basis (ctka,Pka) consists of new fermionic QP 
operators in the E± bands which are hybrids of the d xz - 
and d yz -orbital. In the literature, E- (E + ) is usually 
labeled as the a ((3) band. 

The coherence factors in ^ can be solved to be 



cos 2 k = \ I 



sin 2 6» k = - 1 



4y + £- 



(4) 



Moreover, the single-particle Matsubara Green's func- 
tions in the normal-state are 

S c c (k, lPn ) = - [ dr e^ r ( CkCT (r) C L(0)}o 
Jo 

cos 2 9k sin 2 Ok 



ip n - E- ip n - E+ 

G° dd (KiPn) = - [ dr e lp " r (d kCT (r)dL(0))o 
Jo 

sin 2 Ok cos 2 Ok 



(5) 



ip n - E- ip n - E + 

g° cd (k,ip n ) = g° dc (k,ip n ) = - ( dr e^" r (c k .(r)dL(0))o 

Jo 

sm20k f 1 1 

2 \ipn-E_ ip n -E^ 

Throughout this paper, we shall choose t\ = — 1, t<i = 
1.3, tz = t± = —0.85, and ^ = 1.54, all measured in units 
of |ti|. These represent a good fit to the first-principle 
calculations for the band structures of LaOFeAs^S 



kcr ka 



e x -fj, e 



xy 
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where [i is the chemical potential and (cj^jd^i) are an- 
nihilation (creation) operators that annihilate (create) 
an electron in the d xz - and dj, z -orbital with spin a and 



(1) 



III. RAMAN SCATTERING 

A. Normal state 

Raman scattering intensity is proportional to the imag- 
inary part of the effective density-density correlation 
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function x(q,r) = (T T [p(q, r), p(— q, 0)]) in the q — > 
limit. For the current two-band system, effective density 
operator associated with Raman scattering can be given 
by 

kcr 

+ 7^4 +q >K.(T)], (6) 

where 7^ (7^) is the Raman vertex associated with the 
electrons on d xz (d yz ) orbital. In the Matsubara fre- 
quency space, symmetry-channel dependent irreducible 
Raman response function in the normal state has been 
solved to be 

k 

' /(£_)-/(£+) f(E+)-f(E-) - 
[iui n - (E + -EJ) ioj n - (£_ - E + )\ ' 

where f(E) = l/[exp(i? /ksT) + 1] is the Fermi distribu- 
tion function and 

r k = sin2# k ( 7 £ + 7 £) (8) 

is a weighting factor for the symmetry-dependent Raman 
scattering of a two-band system. T k is comprised of two 
parts where sin2#k corresponds to the effect of orbital 
hybridization, while ( 7 £ + 7k ) is the sum of the Raman 
vertexes associated with two individual orbitals. In the 
case of zero orbital hybridization (e xy — > 0), sin26>k — > 
and hence xn — > 0. 

When the energy of incident light is much smaller than 
the optical band gap of the system, the contribution from 
the resonant channel is negligible. Consequently, Ra- 
man vertex can be obtained in terms of the curvature 
of the band dispersion, known as the inverse effective 
mass approximation^ That is, depending on the sym- 
metry of the Raman modes, ('jQij ~ d 2 e x /dliidkf and 
(7k)y ~ d 2 s y /dkidkf with k; and kf being the wavevec- 
tors of incident and scattered lights respectively. In the 
current two-band model of iron-pnictides, it is obtained 
that 

(7k)si s ~ \{ E T ~ e l V ) =t lC osk x -t 2 cos k y , 
(7k)si s ~ \{^ x ~4 V ) =hcosk x -h cos ky, (9) 
(7k)B 2s = (7k)s 2s ~ £ x x = £ x y v = 4t 3 swk x sinfc y , 

where e% x = d 2 e x /dk x dk x , etc. 

In Fig.[T] we first examine the behaviors of the weight- 
ing factor r k . Respectively in Fig. [Tfa)-(e), we P^ * 
sin20 k , (7£ + 7k)s2 SJ (7k + 7k)si s > r k ,s 2g , and r k ,s lg in 
the first quadrant of the BZ. As shown, both sin2^k and 
(7k + 7k%2 3 are peaked at (tt/2, tt/2), while (7£ + 7k)i?i 9 
is peaked at (tt, 0) and (0, it). Consequently, Fk,B 2B is 
strongly peaked at (tt/2, 7r/2), while due to the mis- 
match of sin 20k and ( 7 £ + 7 k )si , ^k.B lg turns out to 




FIG. 1. (Color online) Panel (a)-(e): To see the proper- 
ties of the channel-dependent weighting factor I\, sin2#k, 
(7k + 7k)s 2g . (7k + 7k)fli fl) Pk,B 2g , and r k , Slg are plotted 
respectively in the first quadrant of the BZ. Panel (f ) is the 
calculated normal-state Raman intensity. 



have two peaks located at (k x ,k y ) = (0.26-7T, 0.74-7r) and 
(0.74-7T, 0.267t). The symmetry mismatch between sin2(9k 
and (7^ + 7 k ) B lg implies that the overall normal-state 
B\g Raman intensity should be weaker than that of the 
Big channel. The reported normal-state Raman intensi- 
ties of Ba(Fei_ a; Co a ;)2As2 seem to be consistent with the 
prediction^ 

Fig- IUf) shows the calculated normal-state Raman 
spectra in both B\ g and Bi g channels. Apart from the 
feature that B\ g Raman intensity is much weaker than 
that of the Bi g channel, £?2 S -channcl Raman intensity 
shows a peak at oj ~ 6.8<i while £>i g -channcl Raman 
intensity shows a peak at ui ~ 1 M\. According to 
Eq. ([7]), the frequencies where the peaks appear are ac- 
tually predictable. When temperature T — > 0, the Ra- 
man peak for each channel can be estimated to be at 
lu = E+(k x ,ky) — E—(k x ,ky) with (k x ,k y ) correspond- 
ing to the maximum of the weighting factor Tk- There- 
fore for i?2 9 -channcl Raman intensity, the peak is pre- 
dicted to be at u) = E + (n/2, 7r/2)-£ , _( 7 r/2, tt/2) = 6.8ti, 
while the -Big-channel Raman peak is predicted to be at 
lo = £ + (0.747r,0.267r) - (0.74tt, 0.26tt) = 7.3ii. The 
above predicted normal-state Raman peaks may provide 
an alternative route to the measurement of the band en- 
ergy scale for iron-pnictide superconductors. 

In iron-pnictide superconductors, the energy scale of 
the nearest hopping t\ is estimated to be about 0.05 — 
0.3eV. Thus the predicted normal-state Raman peaks 
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could have energy to ~ 0.35 — 2. leV. The lower-bound en- 
ergy (0.35 eV) should well be in the non- resonant regime, 
while the upper-bound energy (2.1 eV) could be in the 
resonant regime and poses a question mark for the above 
results to be valid. Current reported normal-state Ra- 
man data on BaFei_ :r Co : i;)2As2 have only be measured 
up to iOQcmT 1 however^ 

B. Superconducting state 

We next study the Raman spectra in the SC state. To 
do so, in a mean-field level one can add a SC pairing 
Hamiltonian: 

H SC =E( A £4X-ki + A^ t /3l k , + h.c.) (10) 

k 

to the diagonalized Hamiltonian Hq in In Eq. (|TU1) . 
the pairing is considered between the long-lived ay_ a QPs 
and between /3k CT QPs only. That is, interband pair- 
ing is neglected. Since decoupled a and j3 bands are 
originated from the coupled d xz and d. yz -orbitals, the 
kind of SC pairing Hamiltonian (fTOj) automatically in- 
cludes both intra- and inter-orbital pairings in the orig- 
inal fermion basis (cko-, ^ko-)' 32 ' 33 The QP excitation en- 
ergy is then given by E kl = (E? + AjJ 2 ) 1 / 2 (/ = a,/?) 
for the two bands respectively. In this section, for conve- 
nience, E- — > E a and E + — > Ep. 

If the pairing originates from the same mechanism, 
most likely ot\ and ai bands will have the same pair- 
ing symmetry. Similarly fi\ and /?2 bands will also likely 
have the same pair ing symmetry. The Raman spectra re- 
ported in Ref. |2(| do not show a clear activation thresh- 
old rather exhibit a finite intensity down to an arbitrar- 
ily small Raman shift. It gives a strong evidence that 
the SC pairing favors an anisotropic nodal gap rather 
than a full isotropic gap such as the s± state or the s++ 
state. Furthermore, the scanning SQUID microscopy 
measurements seemed to exclude the spin-triplet pair- 
ing states and suggested that the order parameter has 
well-developed nodes [l4|. Within the anisotropic and 
nodal scenarios, the possible candidates are the extended 
s-wave and <i-wave states^ 

After a lengthy derivation, the irreducible Raman re- 
sponse function in the SC state is solved to be 

Xs (q->0,r) = - £ ( 7 »') 2 [^(k,T)^(k,-T) 

+ Ji(k,r)7i,(k,-r)], (11) 

where Gi and Ti are the usual normal and anomalous 
Green functions for a superconductor associated with 
band The intra- and interband vertex functions are 
solved to be 

7 £ Q =7kCos 2 # k + 7 £sin 2 k , 

Ik = 7k sin 2 k + 7k cos 2 # k , (12) 

7 f = 7 f* = sin 20 k ( 7 £ + 7^/2, 




FIG. 2. (Color online) Panel (a)&(b): Schematic plot of the 
a- and /3-band Fermi surfaces and the k-dependent ampli- 
tude of (a) extended s-wave gap and (b) d^-wave gap. Panel 
(c)&(d): Momentum dependence of the effective Raman ver- 
tices of (c) (7k a )s 2g and (d) (7k a )si g - Here we only show 
(7r)s 2g and ( 7 r) Blg because ( 7 f )s 29 = (tD^ and 
(tHbi, - (7D Blg (see text). 

where 7 k and 7 k were defined in ([9]) for both B\ g and £>2 g 
channels. Frequency and channel-dependent Raman in- 
tensity is proportional to the imaginary part of the effec- 
tive density-density correlation function (jTTJ) transformed 
to the Matsubara space. As shown in Eq. (|11|) . Raman 
spectra are contributed by both intraband and interband 
transitions. Nevertheless, due to the little nesting effect 
occurring across different bands, interband transitions 
are negligibly small for the Raman intensity. Thus one 
can safely ignore the interband Raman scattering in the 
present case. Consequently, at T — > 0, Raman intensity 
is proportional to 

m= y d£) 2 (^) — , (is) 

J^p \2EuJ ( W -2 J Bk0 2 +r 2 ' 

where T is the broadening which is set to be 0.08 in our 
calculations. 

As is well-known, Raman scattering is a directional 
probe for SC QP excitations. In the present two-band 
iron-pnictide superconductors, the directional selectivity 
is dependent of two factors [see (|T3"])]. One is due to the 
Raman vertex 7k and the other is due to the symmetry of 
the pairing gap A k . The overall Raman will also depend 
on the detailed locations and topology of the Fermi sur- 
faces of the system. It is worth noting that, as shown ex- 
plicitly in (|13[) . Raman intensity is directly proportional 
to the gap maximum of the pairing gap. 

It is important to first observe how the directional se- 
lection occurs for the iron-pnictide superconductors. We 
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FIG. 3. Frame (a) and (b): Comparison of theoretical fitting 
and experimental data of the Bi g Raman spectra. The SC 
gaps tested are extended s-wave Ak = Ao cos k x /2 cos k y /2 in 
(a) and extended d-wave Ak = Ao sin k x sin k y in (b). Frame 
(c) and (d) are the calculated Bi g Raman shift using the same 
parameters in frame (a) and (b). 
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FIG. 4. Plot of the integrand, (7 k i ) 2 (|Ak|/2i?k0 2 , m E q. dHJ) 
in the first BZ. Frame (a)& (b) correspond respectively to I = 
a and /3 with extended s-wave gap; frame (c)& (d) correspond 
respectively to I = a and j3 with d xy -wa,ve gap. 



consider the unfolded BZ for the case of one Fe/cell. The 
cv-band FSs of the 2-orbital model are hole Fermi pockets 
given by i?_(k/) = which are around the T point and 
the corner, (±7r, ±7r). The /3-band FSs are electron Fermi 
pockets given by E + (kf) = which are around the M 
point. In the SC state, the coupling of the two orbitals 
results in complex Raman vertices given in (|12[) . Shown 
in Fig. [5] (c) and (d) are the momentum dependence of 
the Raman vertices (7^") B 2g an d (7k Q )si 9 • Note that for 
the current two-band model, (7^)^29 = (7k Q )s2 S an d 
(7k^)Bi 9 - (7k")si s - :t is because (7k)s 2g = (7k)s2 9 

and (7k)si 9 - (7k)si 9 [ see E q s - © and ®] for the 
current parameters. As shown in Fig.[2](c) and (d), B 2g 
vertex peaks at (±7r/2, ±7t/2) in the unfolded BZ and the 
peak is roughly of same distance to both a- and /3-band 
FSs. As a matter of the fact, it couples roughly equal to 
both a and /3 bands. In contrast, the B\ g vertex has a 
d x 2_ y 2 symmetry and is centered around the T point and 
the corners, (±7r, ±7t). Thus it couples predominantly to 
the a bands in the unfolded BZ. Consequently B\ g Ra- 
man scattering mainly excite the QP in the a bands and 
can give more information about the pairing symmetry 
in the a bands. 

We now test the possible extended s-wave and d-wave 
pairings. Our approach is the following. We will try to 
fit the currently available B 2g Raman intensity^ for each 
possible pairing symmetry. The best fitting parameters 
of the gap magnitudes will be quoted. Using the best 
fitting parameters for B 2g Raman spectra, the predicted 
Bi g Raman intensity will be given. While B\ g Raman 
intensity is also reported in Ref. (2p| . a strong phonon 
mode has appeared at uj = 214 cm -1 responsible for 



Fe vibration, which makes the fitting unfeasible at the 
present time. 

In our following calculations, the only fitting parame- 
ters are A Q and A Q which are in units of t\. Both A Q 
and Ap are adjusted to obtain the best fitting for the 
experimental B 2g spectra. In particular, the B 2g peak 
obtained through fitting is identified with the ratio of 
oj/2Ap which in turn is compared to the actual experi- 
mental data of uj = 69cm~ 1 . One thus obtains the fitting 
value of A^. With the knowledge of A^, one can further 
deduce the fitting values of A Q and t\. The results for 
all possible pairing candidates associated with extended 
s-wave and c?-waves arc listed in Table Q] 

We first consider the case of the extended s-wave pair- 
ing: A^ = A; cos(fca;/2) cos(k y /2) (I = a, (3) shown in 
Fig. [2ja). In a close observation of the B 2g Raman spec- 
tra reported in Refi^, it is identified that only a sin- 
gle peak develops at uj = 69 cm -1 . This is an impor- 
tant point in terms of theoretical fitting. In our fit- 
ting, the key is thus to ensure that both bands give a 
peak at the same frequency (w = 69 cm -1 ). This one- 
peak scenario of fitting may seem unrealistic, but indeed 
it is the only way to successfully describe the presently 
available data. It is found that for the best fitting [see 
Fig. OJa)], A a / Ap = 0.35. The Raman peak occurs 
at uj/2Afj = 0.38 which corresponds to gap magnitudes 
Ap = 91 cm -1 and A a = 32 cm -1 . Moreover we have 
obtained ti = 455cm~ 1 . 

In view of Fig. [3{a) for the best fitted curve, it is seen 
that /3-band contributes most to the overall B 2g Raman 
intensity. As mentioned before, B 2g Raman vertex cou- 
ples roughly equal to the a and f3 bands, thus the indi- 
vidual contribution to the Raman intensity will depend 
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crucially on the multiple effect of Raman vertex, sym- 
metry and magnitude of the gap function, as well as the 
Fermi surface topology. To see this multiple effect, we 
have plotted in Fig. H the integrand, (7^) 2 (|A(J/2.E k i) 2 , 
in Eq. (|T3| in the first BZ. It is shown that /3-band has 
much stronger intensity near the /3-band FS which results 
in a much stronger contribution to the overall Raman 
scattering [see Fig.^a)]. Moreover, due to the nature of 
an extended s-wave gap which has nodes at the BZ edges, 
Raman intensity is linear dependent at small frequencies. 

The case of the d^-wave pairing in the unfolded BZ: 
= A; sin(fc x /2) sin(/c y /2) (I = a,/3) is studied next. 
As shown in Fig. EJb), the gap amplitude in ai FS is 
larger than that in a,\ FS. Thus the Raman intensity due 
to the a-band will mainly contributed by the a2-band 
FS. Due to the same gap amplitude near /3i and /3 2 FSs, 
they will contribute equally to the Raman intensity. In a 
similar approach, to obtain the same Raman peak for the 
two-band model, we again use the two-gap approach. It 
is found that A a /Ap = 0.3 will give the best fitting for 
the experimental Bi g data. Moreover the B^ g Raman 
peak corresponds to w/2A^ = 0.3 which in turn gives 
Ap = 115 cm -1 , A Q = 34 cm -1 , and t\ = 575cm -1 . 
While a band FSs are fully gapped, /3-band FSs are 
gapped with a node. Therefore the Raman intensity is 
powers-law dependent at low frequencies. Moreover, as 
/3-band contributes most to the Raman intensity because 
of larger gap size, the low-energy Bi g Raman intensity is 
actually linear. 

The fitting to the Raman response in B\ g channel 
is not feasible at the moment. There occurs a strong 
phonon mode at 214 cm" 1 due to the Fe vibration^. 
The phonon mode is expected to be removed by chang- 
ing the crystal structure slightly. We have calculated 
and predicted the B\ g Raman responses for the extended 
s-wave pairing symmetry shown in Fig. |3jc) with the 
same parameters as those used in Fig. EJa). As shown 
in Fig. [2£d), although B\ g mode couples predominantly 
to the a band in XM directions, the gap amplitude of a- 
band is smaller than that of the /3-band. Consequently it 
results QP excitation from the two bands with the same 
weight approximately The low-energy Raman response 
is predicted to be power-law dependent due to the full 
gap in both FSs. 

We have also calculated and predicted the B\ g Ra- 
man responses for the d xy -w&vc pairing symmetry shown 
in Fig. |31^d) with the same parameters as those used in 
Fig-CUb). As shown in Fig.^b), cti FS is near gap node 
while a2 FS is fully gapped. Thus ai-band will con- 
tribute more to the low-energy Raman scattering than 
the a2-band. Since the gap amplitudes in /3i FS is equal 
to that of fa FS, they have the same contribution to Ra- 



man scattering. Although the gap amplitude of a-band 
is smaller than that of the /3-band, B\ g mode couples pre- 
dominantly to the a band in XM directions however. It 
results that a bands contribute more to the Raman scat- 
tering than the /3 band. The low-energy Raman response 
is predicted to be power-law dependent with frequency 

TABLE I. Summary of the fitting results for the three pos- 
sible pairing symmetries on the Big Raman intensity.— Row 
1 corresponds to the energy ratio of the Bi g peak (occurs at 
69cm - ) to 2Ay3. Row 2 and 3 are the gap amplitudes of A^ 
and A Q in units of cm -1 and row 4 are their ratios. Row 5 
are the deduced values of t\ in units of cm - . 

Pairing Symmetry extended s-wave d xy -w&ve d x 2_ y 2-wa,ve 



lj/2A 0.38 0.3 0.24 

A^cm- 1 ) 91 115 143 

A^cm- 1 ) 32 34 14 

Ac/A/3 0.3 0.3 0.1 

^(cm -1 ) 455 575 7150 



due to the full gap in a\ band FS. 

To complete the studies, we have also studied the case 
of c? a; 2_ 1) 2-wave pairing symmetry. It is found that the 
best fitting of the B^ g Raman intensity is given by the 
ratio A a /Ap = 0.1 which in turn gives t\ = 7150cm -1 . 
Based on the given unrcalistically large t\, we conclude 
that d x 2_ v 2--wewe pairing is ruled out in terms of the 
current available Raman scattering data. Table U is a 
summary of the the fitting results for the Bi g Raman 
intensity. 

IV. SUMMARY 

In summary, we have studied the Raman response 
of iron-pnictide superconductor in both normal and SC 
states based on a two-band model. Predictions are given 
for the normal-state Raman intensities. A more quan- 
titative fitting to the currently available Raman spectra 
is made to which useful fitting parameters are quoted in 
terms of the gap amplitudes on both bands. 
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